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ABSTRACT 

The interactions between longitudinal vortices and accompanying waves considered here 
are strongly nonlinear, in the sense that the mean-flow profile throughout the boundary 
layer is completely altered from its original undisturbed state. Nonlinear interactions be- 
tween vortex flow and Tollmien-Schlichting waves are addressed first, and some analytical 
and computational properties are described. These include the possibility in the spatial- 
development case of a finite-distance break-up, inducing a singularity in the displacement 
thickness. Second, vortex/ Rayleigh-wave nonlinear interactions are considered for the com- 
pressible boundary-layer, along with certain special cases of interest and some possible so- 
lution properties. Both types, vortex/Tollmien-Schlichting and vortex/Rayleigh, are short- 
scale/long- scale interactions and they have potential applications to many flows at high 
Reynolds numbers. Their strongly nonlinear nature is believed to make them very relevant 
to fully fledged transition to turbulence. 


1 Research was supported in part by the National Aeronautics and Space Administration under NASA 
Contract No. NAS 1-1 8605 while the authors were in residence at the Institute for Computer Applications 
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


The majority of any transition process from laminar to turbulent flow poses a consider- 
able theoretical challenge since fully fledged transition is a strongly nonlinear process, i.e., 
it completely alters the mean-flow profile from its original laminar state. Much interest has 
therefore arisen recently in nonlinear three-dimensional (3D) interactions between longitu- 
dinal vortices and accompanying nonlinear waves, because of the strongly nonlinear nature 
of vortex/wave interactions and hence their likelihood of increased relevance to fully fledged 
transition, e.g., as observed experimentally by Klebanoff et al. (1962), Hama and Nutant 
(1963), Nishioka et al. (1981), Nishioka and Asai (1984, 1985), Kachanov and Levchenko 
(1984), Williams et al. (1984, 1987), Thomas (1987) and simulated computationally by 
Wray and Hussaini (1984), Kleiser and Schumann (1984), Gilbert and Kleiser (1986, 1988), 
Zang and Hussaini (1986, 1987), Fasel et al. (1987), Spalart and Yang (1987), Laurien and 
Kleiser (1989), Zang and Krist (1989). In this work vortex/wave interactions are considered 
theoretically in the context of boundary-layer transition. 

There appear to be three main theories developed so far, of a truly nonlinear nature, 
in the case of flat-surface boundary layers: first, for nonlinear Tollmien-Schlichting (TS) 
interactions via triple-deck-like theory (e.g., Smith (1979, 1985, 1988, 1989), Hall and Smith 
(1984)); second, for nonlinear Euler interactions (Smith and Burggraf (1985), Smith and 
Stewart (1987), Smith, Doorly, and Rothmayer (1989)); and, third, for nonlinear vortex/wave 
interactions (Hall and Smith (1987, 1988, 1989), Bennett, Hall, and Smith (1988), Smith and 
Walton (1988, 1989), Bassom and Hall (1989)). The third, which is our concern here, can be 
regarded as stemming from the first, as Smith and Walton (1988, 1989) show. Alternatively, 
in the presence of curvature, longitudinal vortices are able to exist without any forcing from 
a wave system (Hall (1982), Hall and Lakin (1988), Hall and Seddougui (1989)). Here it is 
found that strongly nonlinear vortex flows completely restructure the boundary layer in a 
manner which renders it neutrally stable to travelling waves; the latter result is closely related 
to the ideas of Malkus (1950) who proposed a ‘marginal theory of turbulence.’ All the above 
theories, for high Reynolds numbers, are known to be inter-connected in that triple-deck 
interactions can lead on to Euler-stage flow (see Smith and Burggraf) and to vortex/wave 
interactions (see Hall and Smith (1987) and Smith and Walton (1988, 1989)), and further 
likely connections are found to arise in this study. All are particularly attractive because they 
can produce substantial alterations of the mean-flow profile or shear in the boundary layer, in 
contrast with linear instability theory and weakly nonlinear theory (e.g., Stuart- Watson- like, 
Benney-Lin, resonant triads, and associated near-neutral analyses) where some interesting 
low-amplitude phenomena can be predicted but nevertheless the mean-flow quantities are 
little changed. Again, the truly nonlinear interactions above are based on rational arguments 
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as distinct from the interesting but ad hoc approaches of certain other theories, while the 
largeness of the Reynolds number taken throughout seems not unlikely to be appropriate 
to the experimental range of greatest concern. The question of which strongly nonlinear 
interaction applies in any particular experimental configuration depends on the amplitudes 
and spectra of the input, at an initial time and/or position. In particular, the present, 
vortex/ wave, interactions apply for waves of small amplitude which nevertheless are able to 
alter the mean-flow quantities because of the relatively short wavelength, compared with the 
relatively large development length of the mean flow, thus inducing full nonlinearity. 

The aim in this work is to describe distinct cases of vortex/wave interaction in boundary- 
layer flows, as well as some interesting and useful sub-cases, of which there are many. We 
summarize the derivations from underlying flow structures, and the governing equations, for 
a number of vortex/ wave interactions, to focus attention on the wide range of possibilities 
and applications, and some solution properties are also given in the form of linear, secondary 
instability, weakly nonlinear, fully nonlinear similarity-type, and nonlinear breakdown phe- 
nomena, along with computational studies, as a start. A main task here however is felt to 
be to put the typical controlling equations of vortex/wave interaction on record, as they are 
believed to be of much significance and of broad application (see below). 

The scales and amplitudes of the vortex flows tend to be relatively simple on the whole 
and can be inferred from the early Taylor-vortex calculations by say Davey (1962), whereas 
the waves* scales and amplitudes are usually more involved. The origins of these scales 
and of the corresponding flow structures are to be found in the flow properties described 
by Hall and Smith (1984, 1987, 1988, 1989), Bennett, Hall, and Smith (1988), Smith and 
Walton (1988, 1989), for lower-amplitude interactions. The two major kinds of nonlinear 
vortex/ wave interaction that emerge, and are addressed below, concern what are effectively 
TS waves and Rayleigh waves. 

Vortex/TS-wave interactions are discussed in Section 2 (see also Figure 1), where the 
scales, flow structure and controlling equations are presented, and in Section 3 where certain 
solution properties are considered. This is mainly for the incompressible boundary layer 
although the corresponding compressible version follows the same pattern, cf. Smith (1987), 
Blackaby (1990), Smith and Walton (1988, 1989). Vortex/Rayleigh- wave interactions (again 
see Figure 1) are then described in Sections 4 and 5 for compressible boundary layers, given 
the earlier groundwork, and given the susceptibility of such boundary layers to Rayleigh 
waves, with the incompressible and other limiting regimes then being obtainable as special 
cases (see Section 5 and Appendix B). A number of interesting nonlinear flow properties 
seem to be suggested. These include the possibilities that the vortex/TS interaction can 
provide a lead into vortex/Rayleigh interaction, even in the incompressible regime, that ei- 
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ther interaction can produce eventually an Euler-stage flow, within a finite distance in the 
spatial problem or within a finite time for the temporal case, that many waves can be acti- 
vated together, and that successive shorter-scale vortex-wave interactions can be provoked, 
producing a cascade of scales. In addition, as the whole boundary layer is changed substan- 
tially in any of these nonlinear interactions, there is potential relevance to by-pass transition 
throughout as well as to more gradual transitions triggered initially by linear disturbances. 
Further comments are presented in Section 6. 

The nonlinear vortex/wave interactions studied here also hold for many other flows in 
principle and have numerous applications. Examples are channel flows (see also Bennett et 
al. (1988)), pipe flows (see also Walton (1990)), wakes, plane Couette flow, water motions, 
e.g., Langmuir circulations, free shear layers, separations, flow over surface roughnesses, 
vortex breakdown, and possibly other rotating fluid flows. The vortex/ wave interaction 
seems to apply in fact to any flow that admits relatively short-scale waves. Among the 
interactions, those of Section 4 for nonlinear Rayleigh-like waves would appear to have the 
broadest application. 

In the following, the large Reynolds number Re is the global one, based on the airfoil 
chord and free-stream speed in the aerodynamic context, as are the corresponding nondimen- 
sional coordinates x, y } z and velocities u } v, w ) streamwise, normal and spanwise respectively 
for flat-surface flow, and the time t . Similarly, the nondimensional density, viscosity, tem- 
perature and pressure are p,/z, T,p (with freestream value poo)> while Moo } a } C denote in 
turn the free-stream Mach number, the Prandtl number, and the constant in the Chapman 
viscosity law which is assumed for definiteness. The characteristic boundary-layer thickness 
is then 0(Re“ 5 ) in terms of p[= Re” 5 p], and the standard TS scalings, i.e., 3D triple-deck 
with y for instance scaled as j, | powers of Re” 1 and x, z as | powers, apply in the setting 
of Section 2, whereas standard Rayleigh scales apply for Section 4. Throughout, the motions 
being considered are 3D and unsteady. 

2. NONLINEAR VORTEX/TOLLMIEN-SCHLICHTING INTERACTIONS 

The vortex/TS nonlinear interactions that we address first here are larger-scale ones, in 
which the mean-flow profile of the entire 0 (Re” 2 ) boundary layer is altered from its original 
laminar form. We start with the wide vortex . The scales involved may be derived from an 
order-of-magnitude reasoning as follows. The small TS waves of typical pressure amplitude II 
say, to be determined, have the triple-deck structure, so that their streamwise and spanwise 
velocity perturbations are of order Re”»/i within the lower deck near the surface, where 
II = Re”* A, and h is assumed to be small. The powers of the Reynolds number present here 
are those characteristic of the triple-deck, and hence of 2D and 3D TS waves, as set out in 
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numerous previous works on TS waves alone. In particular the z scale is comparable with 
Re - *. Forcing of the near-surface vortex flow then occurs at the amplitude-squared level, 
including a vortex spanwise velocity of order Re - *^ 2 . This velocity grows logarithmically at 
the edge of the lower-deck sublayer, as described by Hall and Smith (1988), and so it is little 
different in the main part of the boundary layer, the main deck. There the typical vortex 
dynamics is controlled by the convective- viscous balance acting on the 0(1) scale in x and 
hence, since the z-scale is of order Re - * , the representative w in the vortex has size 0(Re - *), 
given that the streamwise velocity is of order unity. So nonlinear interaction arises when this 
size of spanwise velocity is comparable with that forced by the presence of the TS waves, i.e., 
with the order Re - */i 2 , from above. The critical size is therefore h ~ Re - *, corresponding 
to the pressure amplitude II = Re - * and confirming the relative smallness of the TS waves, 
which contrasts with their substantial impact on the main boundary- layer motion. The same 
estimate of size results from direct extensions of our previous work on related vortex-wave 
interactions (see references). The back-effect, of the induced longitudinal vortex motion on 
the waves themselves, is felt through the streamwise skin friction which helps to control the 
waves’ response within the sublayer. 

The expansions of the flow solutions in the lower, main, and upper decks may then be 
set down, accounting for the extra logarithmic terms necessary in view of the logarithmic 
sublayer behavior of the vortex spanwise velocity mentioned above. The vortex velocity has 
the form 

[u, v, iu] = 0 [1 , Re - *, Re - *] (2.1a) 

across the majority of the boundary layer, while the wave velocity and pressure are 


u,v,w,p] = 0[Re *£], C— ^-foRej 

0[Re - *, Re - *, Re - *, Re - *]£, 
0[Re - *, Re - *, Re - *, Re - *]£, 


(2.16) 

(2.1c) 
(2. Id) 


in the upper, main, and lower decks respectively. See Figure 1. The resulting governing 
equations for the vortex/wave interaction in this context are therefore: 


subject to 


( vortex ) u x + v-g + w z = 0, 

u t + u u x + v Uy + w u z = -p'(x) + % v, 
w t + uw x + vwjj + ww z = 0 + Wy y , 

U = v = 0, w - -=|(|P| 2 + a~ 2 \P z \ 2 ) at y = 0, 
A 


(2.2a) 

(2.26) 

(2.2c) 


(2. 2d) 
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u — * u e (x),vy —* 0 as y —> oo ; (2.2e) 

( wave ) P Z z ~ X-=-Pz - o?P — QA, (2.3a) 

A 

with 

T = Uy at i/ = 0. (2.36) 

Here the scaled vortex velocity ( u,u,uJ ) and skin friction A, and the wave pressure P and 
displacement —A f are all unknown, as is the (real) wavenumber a(x) if the spatial devel- 
opment is under consideration. Also, p'(x) = ~u e it f e (x) is the prescribed external pressure 
gradient, if any, and the spanwise scale has z = Re“^Z. The wave pressure-displacement 
(P-A) law required to complete the system is given by solving 

(dj + d 1 2 3 z -a 2 )p = 0, (2.4a) 


with p — * 0 in the farfield, (2-46) 

p —> P, Py — > -a 2 A, as y — > 0+, (2.4c) 

for the upper-deck response in the potential flow just outside the boundary layer, in the 
present incompressible regime; and the coefficients appearing in (2.3a) are defined by 


3 A , \ 

2 + 2Ai(£ 0 ) { + Ai'(( 0 )J 


,g = (iaX)lAi'U 0 )/K, 


1 — 7 f°° 

£ 0 = — isfl/(aA)», K= (2.5 a — d) 

where Ai is the Airy function and, for the spatial case again, Q (real) is the constant imposed 
frequency. 

The wave-forcing of the vortex motion appears in the effective spanwise surface velocity 
in (2. 2d), an amplitude-squared effect anticipated in the first paragraph of this section. The 
vortex-forcing of the wave, on the other hand, is through the skin-friction factor ~X in the TS 
pressure equation (2.3a). In the pure spatial- development problem the terms Ut,w t in (2.2b,c) 

3 

are to be omitted, and, in the background, multiple scaling d x — » Re«£?j>r + d x operates, with 
the short-scale wave being dependent on X through the form P( x y Z) exp (iaX in effect 
while the long-scale vortex flow is independent of X [more precisely, iaX should be replaced 
by i J adxRe®]. The main alternative of pure temporal development has d x identically zero 
instead, in (2.2a-c), along with the multi-scaling d t — * Re*df + d t where the fast-scale wave 
varies with t as above, i.e., we have P(t, Z ) exp(zaX — ifR) in effect, but the slow-scale vortex 
evolution is independent of f , and a is then constant with f l(t) to be determined. 

The wide-vortex/TS interaction, then, is governed by the equations (2.2) - (2.4). 
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The next two types of vortex/ wave nonlinear interaction addressed are for the square 
vortex T for which the y- and z-scales of the vortex are comparable, of order Re“*. These 
types have scales implied to some extent by Bennett and Hall (1987) and by an examination 
of the wide- vortex case for condensed \z\ scales. First, the square vortex has the scales 


[ti, v y w] = 0[l,Re 5 } Re *], p - p(x) + 0(Re *), (2.6a) 

with z = Re“*z and H of 0(1), and the wave disturbance is now condensed inside the 
boundary layer, such that its velocity and pressure fields have 

\u,v,w,p) = 0[Re - 5]£i,£i = Q/n Re^ , (2.66) 

0[Re-J, Re"?, Re"?, Re-5]£ a (2.6c) 

2 

in the midst of the boundary layer and in the 0(Re“5) sublayer, in turn. So here the vortex 
equations are 

u x + Vy + w Y = 0, (2.7 a) 

U t + U U x + V Uy + W Uj = —p'(x) +Uyy+Ujj, (2.76) 

v t + u v x + V Vy + w Vj = -P 2 y +Vyy + Vjj, (2.7 c) 

w t + u w x + 57 Wy + w w T = -p 2 z +Wyy + Wjj. (2.7 d) 


These are the nonlinear Gortler- vortex equations in a boundary layer at zero Gortler number 
(Hall (1988)); in the presence of curvature a term proportional to u 2 must be inserted in 
(2.7c). The boundary conditions appropriate to (2.7a-d) are 

u — v — 0,uJ = — =|[|P | 2 + o.~ 2 \P^\] at i/ = 0, (2.7e) 

A 

u — > u e (x ) as y — * oo, (2-7/) 

while the coupled wave-pressure equation is 

Pjz-^jPj-ct 2 P = 0. ( 2 . 8 ) 

Here X is the unknown skin friction as defined in (2.3b) and 7F is as in (2.5), but now the 
interaction relation of (2.4) is replaced by A —* 0 effectively. This new feature is due to the 
shortened streamwise and spanwise length scales associated with the wave, both of which are 
now 0(Re - 5), thus suppressing the inviscid pressure feed-back from outside the boundary 
layer. The other new features of the present two-tier structure, for the square- vortex case, are 
the balance of normal and spanwise diffusion, in (2.7b-d), and the interpretation of (2.7c,d) 
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as a streamwise vorticity equation (on elimination of P 2 by cross- differentiation). It should 
be noted however that (2.8) has no solutions with a real if Ay = 0, and even with Ay = 0 the 
possibility exists that there are no solutions of the full interactive equations in this case. 

The other type of nonlinear interaction involving a square vortex has some analogies with 
Bennett et al.’s (1989) for channel flow. Here the vortex exhibits scales as in (2.6a) again, 
but the wave form is distinct from (2.6b,c) in that the triple-deck structure is reinstated, 
such that 

[u, v,w,p) = 0(Re~«)£ + 0(Re -y )£ (2.9a) 

in the upper deck, and so on; in particular 

p= (Re-iP 0 + Re~*P 1 )£ + ... (2.9 b) 

with Pq(x ) uniform in ~z but Pi depends on x,~z. The relative contributions of order Re - * 
in the wave form affect both the vortex and the TS response throughout. Thus the vortex 
equations are again (2.7a-f) except that the spanwise slip condition in (2.7e) is replaced by 

w = + C.c. + a- 2 |i\y| 2 ] at y = 0, (2.10) 

A 

and the wave-pressure equation becomes now 

Pn y - phpu - a 2 P 0 = GAo (2.11) 

A 

where Aq = a~ 1 Po. Here again the y- and z- diffusion effects in the vortex motion are 
comparable, but the increase in the skin-friction variation across the span, relative to the 
wide-vortex case, causes the splitting of the TS wave response as in (2.9b). The solution of 
(2.11) for Pi can be written in integral form as in Bennett et al. and substituted into the 
slip condition (2.10) in principle. Other scales and forms of splitting are also possible for 
this square-vortex type. 

We move on now to the fourth type of nonlinear interaction to be considered, that for 
a small vortex . This arises as an interesting sub-case of (2.2) - (2.4)’s interaction, with the 
characteristic y-scale reduced to 0(A) say, near the surface, where the new parameter A < 1. 
So we might expect that u ~ A, in an 0(A) sublayer, with the x-scale reduced to the order A 3 
to preserve the convective- viscous balance in (2.2). Hence by continuity 17 is large, 0(A -1 ), 
and w is still larger, 0(A~ 2 ), provided the Z-scale remains intact. The nonlinear interplay 
between the vortex and the wave also stays intact if the scaled wave pressure P is increased 
by an amount A -1 , from (2.3), with a remaining 0(1). So here the vortex quantities are 

[u,v, w\ = 0[A,Re - ^ A -1 , Re - * A -2 ] with jF ~ A, x ~ A 3 , (2.12a) 
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and the TS wave pressure is increased to 

p = 0(Re - 8 A -1 ), (2.126) 

to within a logarithmic factor. The governing equations for this case, in scaled form, are 
(2.2) with p'(x) absent and with (2.2e) replaced by 

Uy — * A a as y — ► oo, (2.13) 

coupled with (2.3), (2.4) again. Here the constant Ai is the skin-friction factor for the undis- 
turbed incident boundary layer. As A tends back towards 0(1) the wide- vortex case of (2.2) 
- (2.4) is approached, whereas for A reduced towards 0(Re _ ») the small-scale interaction of 
Smith and Walton (1988, 1989) is recovered. A number of other small- vortex/ wave nonlinear 
interactions can be derived in similar vein from the previous ones. 

All the vortex/wave interactions described above are fully nonlinear in the sense defined 
in Section 1. In general each one needs a computational treatment, marching forward in x 
from given starting conditions at x = 0, say, in the spatial-development setting, or forward 
in t from initial conditions at t = 0 for the temporal setting. Note however that for either 
the temporal or spatial problem we cannot arbitrarily choose the input vortex velocity field 
since the associated Tollmien-Schlichting wave must satisfy an eigenrelation dependent on 
the initial shear stress. Some solution properties are presented below. 

3. NONLINEAR INTERACTION PROPERTIES 

Most of our interest here is in the spatial development for the wide vortex of (2.2) - 
(2.4) and for its small-vortex form (2.13). The nonlinear flow properties of the vortex/wave 
interactions seem to depend to a large extent on the amplitude and spectra of the input 
disturbance upstream. 

Certain special cases may be addressed first, as guidelines. Thus if the input comprises 
two oblique waves of relatively low amplitude and spanwise wavenumbers i/3 say, then an 
analysis for (2.13) for instance can be conducted similar to the analysis in Hall and Smith 
(1988) (but corrected for a logarithmic effect, like that in (2.1b), as described in a subse- 
quent paper with Professor P. J. Blennerhassett). Such analysis yields interaction equations 
between the near-neutral waves’ amplitudes and the induced vortex motion. The solutions 
of these weakly nonlinear equations in the Hall and Smith (1988) case show that, depending 
on /3, either a finite-distance singularity is encountered or a far-downstream asymptote with 
exponential growth is attained. In the present context, the former would tend to reinstate 
the full vortex/wave system associated with (2.13), whereas the latter would lead on even- 
tually to the longer-scale vortex/wave interaction of (2.2) - (2.4). Similarly, if the input 
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wave is near-planar then analysis along the lines of Smith and Walton (1988, 1989) can be 
applied initially. For relatively small amplitudes, the solution properties found with weak 
nonlinearity present again tend to reactivate (2.2) - (2.4) or (2.13) in full. For non-small 
amplitudes, pronounced secondary instability to 3D modes is found, among other things, 
with the near-planar input, the 3D components growing initially in a form analogous to the 
exponential of an exponential of distance. 

The above cases could be used to provide upstream starting conditions for the full 
vortex/ wave- interact ion systems (2.2) - (2.4), (2.13), e.g., at finite or large negative x re- 
spectively. The ultimate downstream behavior on the other hand for the small-vortex case 
is considered in Smith and Walton (1988, 1989), where three main possibilities are raised. 
One is that the full interaction continues to downstream infinity and acquires a nonlinear 
similarity form (with the scaled u, P, y, W behaving as scaled x» , x - ?, s?, x“j in turn), prop- 
erties of which are given in the last-named paper and by Walton (1990). In that event, 
the small- vortex/ wave interaction acts as a precursor to the wide- vortex form. The second 
possibility concerns a 3D-strong-attachment singularity occurring in the vortex and wave 
solutions at finite x downstream. This typically takes the form, with n > 

V ~ { x a ~ x) n V,v ~ ( x A - x)~ n ,w ~ {x A - x)~ 2n , (3.1a) 

on approach to a surface attachment line x — x A (Z)—, say. Here u is typically 0(1). Near 
a valley plane Z = Z\ for instance, where w ~ {Z — Zi)w(y) vanishes, the cross-stream 
balances 

v= + w = 0, vw= + w 2 = w= = 

dominate, yielding 

v = 7[exp(-7y) - l],ur = 7 2 exp(-7f), (3.16) 

where 7 is a positive constant. Hence, the wave-pressure amplitude P is proportional to 
(x A — x)~ n and also becomes singular at x = x A . The third possibility in Smith and Walton 
is again a finite- distance singularity but of a 3D-separation kind, taking place in the vortex 
motion, probably at a peak plane Z = Z 2 say. There 

V ~ { x b — tZ ~ ( xg — x) q U, (3.2a, 6) 

v ~ ( xb - X )~ T , w ~ ( xb - x) q ~ 1 (Z - Z 2 ), (3.2c, d) 

with N = T + q — 1 positive and q > 1 — 2I 1 , so that the local response is predominantly 
inviscid in a region which thickens in singular fashion as x — * Xq — - The solution has the 
form 

dU/dY = AiU mi + A 2 U Tn7 y m\ = 1 + N/q,m 2 = mi — l/q (3.2e) 
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where Ai,A 2 are constants. For example, if T = 5, q = 1, then N = | and U takes a (tan) 2 
form. The description (3.2a-e) holds for a finite range of Y, with U becoming singular at 
Y = Yi say in such a way that a thinner region of extent 0(1), in terms of y — ( xb — x) _Ar Yi, 
is induced to smooth out the velocity profiles. In this 0(1) region, the velocities u, w are 0(1) 
and their profiles are arbitrary, i.e., dependent on the flow history, apart from the matching 
with (3.2e) at the lower extremes and (2.13) for example in the upper extremes. In particular 
the effective boundary-layer displacement 

6* ~ (x B - x)- N (3.2/) 

becomes large in this separation singularity, connected perhaps with the formation of a 
lambda (loop) vortex in practice. 

The proposed attachment singularity (3.1) and separation one (3.2) apply equally well to 
the wide- vortex case in principle. The similarity form mentioned previously does not apply, 
however, and appears at first sight to be replaced by the far-downstream behavior 

P ~ P 0 x~* + P\X~% + . . . ,a ~ a 0 x~*, (3.3a, b ) 

u ~ 1,17 ~ x“s ,W ~ sc -1 , A ~ x - *, (3.3c — /) 

suggested by (2.2) - (2.4), at least if the Z-scale stays fixed at 0(1), e.g., with Z-periodicity 
present. Here Pq is independent of Z , and the splitting in (3.3a) and the rest of the solution 
structure bear resemblance to the square- vortex form in (2.9)-(2.11), although now |£ 0 | is 
large (~ ij), so that relatively high-frequency features apply. Conversely, a low-frequency 
input with fi < 1 would produce the form in (3.3), downstream, with x replaced by fi -2 
in effect. Further investigation however suggests that there are no high-frequency solutions 
with \£o\ large, e.g., from analysis of (2.3a) with (2.5a-d). Instead, the flow solution seems 
more likely to terminate with a finite-distance singularity, as in (3.1), (3.2) or with the wave- 
pressure amplitude tending to zero in a square- root fashion, in view of (2. 2d), similar to a 
special case given in Smith and Walton (1988, 1989). Other possibilities for the wide-vortex 
and small- vortex nonlinear interactions with TS waves may exist of course, and likewise for 
the square-vortex forms in Section 2. 

Computational studies have been made of both the wide- and the small- vortex/ wave 
interactions of (2.2) - (2.4) and (2.13) respectively. Spectral treatments for the Z- variation 
have been applied to each case, and a finite- difference Z-representation has also been applied 
to the small- vortex case. These studies are being continued. Here we report on some prelim- 
inary but significant results obtained in our numerical investigation of the wide-vortex/wave 
interaction problem. This was done using a mixed finite-difference/spectral approximation 
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to the vortex/wave equations (2.2) - (2.4). Thus, for example u is written as 

oo 

u = Ti 0 (x,y) + ]T)u n (x,$r)cos n/3Z, (3.4) 

l 

and the x,y dependences of u n are then approximated using finite differences. All of our 
calculations were for the zero pressure gradient case, p' = 0, but the scheme as described 
could be carried out for pressure-gradient-driven boundary layers. 

Suppose then that u,v and w together with P and A, and a,Q, are known at x = x; we 
now describe a scheme which can be used to advance the solution to x = x +? in such a way 
that the wave frequency Q, is held fixed. Thus we are assuming that as the wave evolves its 
frequency stays constant whilst its wavelength and amplitude vary. 

In order to step the solution forward, we decouple the vortex and wave equations by 
making w at x = x + 1 satisfy the required boundary condition at y = 0 evaluated in terms 
of P known at x = x. The x and z momentum equations for the combined mean flow- vortex 
field can then be stepped forward using essentially the scheme used by Hall (1988) in an 
investigation of fully nonlinear Gortler vortices. The reader is referred to that paper for 
precise details of that scheme. It suffices here to say that y derivatives are approximated 
using central differences and the nonlinear terms involving any harmonic content of u,v and 
w are iterated upon until a converged solution is found. 

The above procedure is used to advance the vortex velocity field to x = x + I and 
the corresponding shear stress at that point can then be found from u. The eigenrelation 
specified by (2.3a) and (2.4a,b,c) will, for fixed frequency, determine a complex value for the 
streamwise wavenumber a. However, it is implicit in our analysis that the wavenumber a is 
real, so that the eigenrelation at x = x + e in general does not have an acceptable solution. 
At this stage there are two simple procedures which can be used to remedy the situation. 
Firstly, the value of w(y = 0) at x = x + e can be iterated upon in order to make a calculated 
at x + J purely real. In effect this is most easily done by iterating upon some measure of 
the Tollmien-Schlichting wave amplitude at x = x + e. Alternatively, we proceed by writing 
(2.3a) in the form 

where — and + denote quantities evaluated at x and x + e respectively. The right-hand side 
of (3.5) is then known, so that after expanding P in a Fourier series (3.5) can be solved for 
P at x + e. The value of a used to determine P is then iterated upon in order to make (say) 
(P)i = 0 at Z — 0. In this procedure the Tollmien-Schlichting frequency is of course fixed 
and a value for a at x + e is obtained. 
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The calculations which we report on have been carried out using the second of the proce- 
dures described above in order to keep a real. The input used to begin our calculations was 
found in the following way. Firstly, we assume some form for u at an initial value x = x* 
and then solve the eigenvalue problem determined by (2.3a), (2.4a, b,c) for the appropriate 
real values of a and ft. The values of a and ft calculated in this way depend on u through 
the wall shear A, so that for example an almost 2D input Tollmien-Schlichting wave can be 
constructed by choosing u n in (3.4) for n > 1 to be small compared with u 0 . Having solved 
for the Tollmien-Schlichting eigenrelation at x = x* the boundary condition to be satisfied by 
w then can be determined in terms of P. The initial profile for w at x = x* is then chosen to 
be consistent with this condition. In the calculations reported here (3 = 0.02, x* = 55, and 8 
Fourier modes were retained in the Fourier expansion of the vortex and Tollmien-Schlichting 
fields. The initial distribution for u was taken to be 


u 



cos n(3Z 



(3.6a) 


where ub is the Blasius profile and, Uq = e v — e 2y and A a parameter which can be varied 
so as to alter the size of the incoming vortex. The spanwise velocity w was then taken to be 

OO 

w = SXnWo(y) sin nfiz (3.65) 

i 

with Wq = cos ye -y and % n chosen so that w satisfies the boundary condition on w at y = 0 
determined in terms of the Tollmien-Schlichting pressure. Clearly the above initial conditions 
are rather arbitrary but this is always the case with longitudinal vortex calculations unless the 
receptivity problem is discussed; see Hall (1983, 1989). Some limited experimentation with 
other initial conditions produced qualitatively similar results but we do not claim to have 
made an exhaustive investigation of the effects of the initial conditions on the vortex/wave 

interactions. . - 

In order to monitor the evolution of the vortex and the wave, the following quantities 
were calculated as the flow was allowed to develop: 

e n=/ (K + w 2 n )dy, n= 1,2,... (3.7a) 

JO 

and P n where 

OO 

P = Po + X] Pn cosn(3Z. (3.7 b) 

i 

In Figure 2 we show results from a calculation where the initial wave is almost two-dimensional. 
Figures 2a, b show the development of e x , e 2 , e 3 , e 4 and Po, Pi, . . . , P 4 respectively. Figure 2c 
shows the corresponding development of a. The calculations were started from x = 55 with 
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more three-dimensional as it moves downstream and that at a finite value x a singularity 
appears to arise, beyond which the solution cannot be calculated. The apparent singularity 
occurred in the same place when the finite- difference resolution was increased or the spectral 
resolution decreased. It was not possible for us to perform the calculations with more than 
eight Fourier modes because of our limited computing resources. However, we believe that 
Figure 2 represents a significant calculation because it demonstrates the three-dimensional 
secondary instability, and ultimate breakdown, of an initially two-dimensional wave and the 
induced vorticity structure. The resolution of the calculations was not great enough to deter- 
mine whether the singularity encountered numerically is related to the ones discussed earlier 
in this section. 

In Figure 3 we show the results of a similar calculation but with smaller initial strengths 
for the vortex and wave fields. Here the vortex and waves decay as x increases. We note 
that the wave appears to effectively disappear before the vortex field. This is consistent 
with there being a subcritical bifurcation of a three-dimensional wave from a longitudinal 
vortex velocity field at a finite value of x (see our earlier comment concerning a square-root 
wave-pressure behavior). 

4. NONLINEAR VORTEX/RAYLEIGH INTERACTIONS 

Nonlinear interactions between longitudinal vortex flow and inviscid Rayleigh waves are 
considered here for the compressible or incompressible boundary layer, with, as a result, the 
boundary layer’s 0(1) mean-flow profile again being changed completely from its original 
form. 

The scales involved may be deduced mostly from a first-principles argument. Thus if the 
induced Rayleigh wave has pressure amplitude 7r, its typical velocity amplitudes are also of 
order W, by its inviscid nature, and the representative wavelengths are all of the short-scale 
size Re - ^. Hence the nonlinear inertial effect provoking a mean-flow correction is of order 
7f 2 Re*, e.g., from uu x , uv x , uw x . This is to be compared with the minimum inertial force in 
the typical long-scale vortex motion, namely Re - *, for a full “square” vortex of size as in 
(2.6a); this force is from the spanwise and normal momentum of the vortex, e.g., uv x ,uw x 
(and the viscous forces such as Re -1 ^ y), rather than the strong streamwise momentum force 
of order unity. So the wave affects the mean flow at zeroth order if ¥ 2 Re 5 is comparable 
with Re - * , i.e., if the wave pressure has amplitude tt ~ Re - * (smaller than in Section 2), in 
principle. There is a complication, however, similar to that in the vortex/TS case, namely 
the appearance of logarithmic behavior in the induced vortex velocities close to the (linear) 
critical layer(s), situated at y = f say. The logarithmic response arises because the 3D 
wave velocities there grow like ( y — /) -1 (see below), so that the nonlinear inertial spanwise 
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wave velocities there grow like (y — /) -1 (see below), so that the nonlinear inertial spanwise 
forcing is proportional to (y — /)~ 2 , which can be balanced only by the viscous term Wyy of 
the vortex if w is proportional to 2n\y — f |. Hence logarithmic contributions are drawn into 
play again, slightly altering the interactive balances. 

In consequence the flow solution in the nonlinear vortex/ compressible- Rayleigh interac- 
tion (see also Figure 1) has the underlying form 


[u,v,w } p,p,T,(i] = [u, Re Ju,Re~*u;,p(x) + Re 1 p 2 ,p,T,p] 
+Re - 5[u^\ v^\ p^\ p^\ T&, + . . . 


(4.1) 


where the first square brackets describe the vortex motion and the second the wave, while 
Ci — (|£nRe)“5 is small. The expansion (4.1) applies across most of the boundary layer, 
with y,z ~ 1, and leads to the compressible vortex equations 


(? u)x + (p u)y + (p w)x = 0, (4.2a) 

p(u U x + VUy+WUj) = - p'(x ) -(- (/I Uy)y + (/Z U*)*, (4.26) 

p(u % + v Vy + w Vj) = —p 2y + (p Uy) x + (2p Vy)y + (/l(Wy + vj))z, (4.2c) 

(4. 2d) 
(4.2e, /) 
( 4-2 g) 

( 4.2 h) 
( 4-2 0 


(4.3a) 

p — » 0 (or outgoing waves) as y — * oo,py = 0 at y = 0. (4.36, c) 

Here the wave has pM = pexp(iaX - ittt) + c.c., with a,Q(= ac ) real, the amplitude p 
is independent of the faster scales X = Re5x,i = Re = Z, and similarly for etc., 

leading to the inviscid response in (4.3). Also, p 1 = —u e u' e is the external pressure gradient, 
and p 2 = PihMlc — (p' — | p)(u x + % + w T ). The nonlinear interaction occurs through the 
definition of M required for (4.3), 

M = (u — c)p 5 Moo, (4.4a) 


p(u w x + vWy + ww z ) = -Pzz + (p u z ) x + {p(wy + Vj))y + ( 2p Wj)j, 

pT = 1, p = CT, 

p(u T x + V Ty + w Tj) - cr~ 1 {{jL Ty)y + (/I T j)j} , 


with 


u = v = w = 0, T = T w at y = 0, 
u — > u e (x), T — > 1, Vy —* 0 as % — > oo; 

and to the compressible Rayleigh equation for the effective wave pressure p , 


Py y Pz Z yPy 4" M z pJ ) (X (1 M )p 0 


with 
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I ii| nil N> r ,v i> , ii i If *■ 


and through the boundary conditions on the vortex flow at the critical layer 

7- 2 /!- 1 A "* 


? = ‘ k. . {« 2 ^(l?| ! ) + A-'^lfel 2 ) - 4 A-VtA tIAI’} (4.46) 

) 2 a 2 


at 


y = f(x, z ) ± [where u = c]. (4.4c) 

Here q is the cross-flow velocity tangential to the critical- layer curve y = f(x,z), so that 

q = A’TZJ, v = fjw, at y = f, (4.4d) 


and A = 1 + /J. The contribution (4.4a) describes the vortex-forcing of the wave, in essence, 
while the cross-flow slip velocity (4.4b) represents the main back-effect of the wave on the 
vortex motion, thus producing nonlinear interaction in which both the mean flow and the 
wave are unknown. 

The main details behind the slip condition (4.4b) are presented in Appendix A. 

The nonlinear vortex/ Rayleigh interaction is given by (4.2) - (4.4), for the compressible 
boundary layer. The pure spatial-development case is shown above, that for pure temporal 
development having u u x replaced by u t , etc., and a fixed surface temperature is assumed 
as an example. Some of the comments in Section 2 apply here as well, and further points 
about the flow structure are the following. First, the induced slip velocity in (4.4b) clearly 
has a connection with that for TS waves in Section 2 corresponding to / = 0 and the high- 
frequency regime (see also Section 6). Second, the relatively thin critical layer surrounding 
y = f is of the linear viscous kind, because the wave pressure amplitude n f is only 0(Re~’ £i) 
(see references), and the thickness is of order Re“ 5 . From (4.3a), (4.4a), the wave response 
nearby as y —* /± is of the form 

p ~ po + spi + s 2 p 2 + s 3 p 3 + . . . (4.5a) 


where s = y — /, in general logarithmic terms must be absent to allow both a, c to remain 
real, and 

~M ~ sMx + s 2 M 2 + . • • , (4.56) 

with the coefficients above being ^-dependent and related by 

piA = poz/z, (4.6a) 

— 2p 2 A - /i 7 pi + 2 vy -- (/«? i - Poz) + Poj z - a 2 p 0 = 0, (4.66) 

M. i 
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(4.6c) 


-2 p 2 jfr - 2p2(fj j - 3A^ + 2^/x) + P!** + ^-(M 2 / 7 - M a )pa 


2M, 


m, 


+p 1 (-a 2 + -^-f T - ,) + [M 2 (poz i - a 2 po) - 2M 2J poj]/M 1 = 0. 

(In principle these give three equations fixing the three terms PojPijib)- Hence it follows 
that the 3D wave velocities are given locally by 


U ™ F-\S 1 + Fq + . . . , 

v ~ G_\S 1 + Go + * * * , 
w ***> H_is 1 + Hq 


(4.7a) 

(4.74) 
(4.7c) 

as anticipated at the start of this section; in addition p ~ — Pzpaz/i'yct 2 M^As 2 ). Here the 
inviscid disturbance equations show that 

l M LPo G -i = -Pi/ioiui, 


'iMIcPqH-i = (fzpi - p^/iau!, 
7 M^p 0 F_! = -H^j/ia, 


(4.8a — c) 


consistent with (4.6), where p ~ Po + Pi s + • . • ,u — c ~ uis + u 2 s 2 + — See also Appendix B 

for the incompressible regime. The singular behavior (4.7) is responsible for the generation, 

at the amplitude-squared level, of the logarithmic-flow effect which leads to the effective slip 

2 

in (4.4b), via the 0(Re - ’) critical layer (see Appendix A). Third, the 3D-wave condition 
in (4.6) can be shown to agree with the generalized inflexion-point condition for 2D or 
3D simple waves on a parallel flow (and with other simpler cases). The latter waves can 
act as triggering mechanisms for the present nonlinear interaction. Fourth, the continuity 
properties of the total mean flow u,u,u7 across the critical layer are worth noting. These 
are that u,u l , vf',v,v',w are all continuous, with discontinuities appearing first in u" 1 , v", w' 
(the prime denotes the normal derivative), and requiring higher-order smoothing within the 
critical layer, which is addressed in Appendix A. Fifth, the wave’s inertial effects on the 
mean flow are felt solely in the slip condition (4.4b) to leading order, the effects in the rest of 
the boundary layer being negligible (just!) due to the logarithmic response described earlier. 


5. SPECIAL CASES AND LIMIT SOLUTIONS 

The vortex/compressible-Rayleigh-wave nonlinear interaction set up in Section 4 poses a 
computational task, in general, which seems a particularly severe one in view of the unknown 
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moving boundary present (in (4.4)) and the coupled partial-differential systems (4.2), (4.3). 
No full solutions have been obtained yet, and appropriate computational treatments are 
only just being considered. For that reason we turn briefly here to certain special or limiting 
properties in (i) - (vii) below, to help provide some possible guidelines and suggestions. 

(i) The linearized version, where w and d t are identically zero in (4.2) and u,v take 
on the compressible Blasius form, say, applies at small wave-amplitudes \p\ and leaves the 
linear compressible-Rayleigh equation (4.3) controlling instability, with M known in advance. 
Solution properties at various Mach numbers are given by Mack (1975, 1984) and Malik 
(1982, 1987) principally. The main result in the present context is that neutral modes exist 
at all Mach numbers (see, e.g., Mack’s (1984) figures) and so can act as triggers for the 

nonlinear interaction of Section 4. See also (iv) below. 

(ii) Weak nonlinearity can also be handled analytically in principle, an example being 
for two input oblique waves of lowish amplitudes. This is analogous to the Hall and Smith 
(1988) oblique-wave/vortex analysis in the TS case, and includes secondary 3D instability 

at the start, as does the full system of Section 4. 

(iii) Wide-vortex/wave interactions similar to those in Section 2 can arise as limit cases 

of (4.2) - (4.4) for enlarged spanwise scales. 

(iv) Special ranges of the Mach number M«, are of theoretical and practical interest, 
including zero M m (see also Appendix B), small Moo, the transonic range M x -> 1, and 
the hypersonic range of large M^. These have connections, in turn, with the vast literature 
on linear incompressible Rayleigh modes in boundary layers, with Gajjar’s (1989) linear 
and nonlinear critical-layer work, with Bowles’ (1989) linear and nonlinear instability work, 
and with recent studies of linear hypersonic-flow instabilities. Concerning the hypersonic 
range in particular, the undisturbed steady 2D boundary layer with no imposed pressure 
gradient itself acquires a two-layered form at large Mach numbers (Bush (1966), Lee and 
Cheng (1969), Stewartson (1964)), with a relatively wide high-temperature layer, wherem y is 
0 (M^), at the upper edge of which is a relatively thin high-vorticity layer with y - M^F(x) 
of order (InMl)*. Here y = M^F(x) denotes the scaled boundary-layer displacement. 
In line with this, the linear instability modes split into two types, the so-called vorticity 
mode concentrated within the high vorticity layer and having the maximum growth rate 
(Brown and Smith (1989)), and the so-called acoustic modes which spread normally across 
the 0 (M£J layer and have smaller growth rates (Hall and Cowley (1989), Brown and Smith 
(1989)). These linear features agree quite well with Mack’s results at large M ^ and suggest 
the two most likely structures of nonlinear vortex/compressible-Rayleigh interaction at large 

Moo, as follows. 
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First, the high-vorticity form of nonlinear interaction occurs where 


y - Ml 0 F{x) + (2r) 1 y [with T exp(r 2 ) = M^J (5.1a) 

and the scalings operating are of the form 


n=l-^V = MlF,* + 2TV,*=-M^W, 

(5.14) 

m 

z = M 00 ( 2T)~h, 

(5.1c) 


~ 1 , T 

(5. Id) 

\ 

These scalings are implied by the Prandtl shift in y, in (5.1a), which 

introduces in effect 

- = 

a large Gortler number F xx (for the high-vorticity-layer flow) that 

is negative for flat- 

= 

surface flow where F oc -\-x? . Thus the vortex becomes relatively K wide. J 

’ It is also quasi-2D 


in the cross-flow plane, since its governing equations become 


(pVh + ( pW % = 0, 

(5.2a) 


P^xx ~ 

(5.26) 

- 

-p{V W$ + W Wi) = -7T, + (71 w 

(5.2c) 


-p{VT, + WT i ) = a-\ 1 iT i) )„ 

(5. 2d) 

= 

pT = l,]i=CT, 

( 5-2c,/) 


p(V + W uu) = (7Zu lfi ) fi , 

(5.2 ? ) 


from (4.2) with (5.1). The coupled wave-amplitude equation is now 



- A „ 2a 

po a Pr“P - 0, 

M 

(5.3a) 


where 



M = (ci - ui)pJ, 

(5.36) 


however, upon suitable scaling in (4.3), while the back reaction of the wave on the vortex 


motion is felt through the slip condition 



w 3C 7 2(u;(j? c )) 2^ (IpI ) 

(5.4) 



at the unknown critical-level curve, from (4.4). Hence, despite the quasi-planar balances 
in (5.2a-f), the streamwise momentum balance (5.2g) in the vortex still exerts influence on 
the nonlinear interaction due to the appearance of in (5.3b) and (5.4). We observe that 
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an alternative type of nonlinear interaction can arise either with sufficient surface curvature 
present, if the resulting Gortler number is equal to M^Fxx + 0(r _1 ), or with negligible 
curvature F xx , since then the vortex can remain “square”, more like those in Section 4 
and below. Again, for small wave amplitudes linearized properties hold, with (5.2), (5.4) 
reproducing the basic boundary-layer solution via the Howarth-Dorodnitsyn transformation, 
while (5.3) then gives the equation addressed by Brown and Smith for which the neutral 
wavenumber is a = 1. Comments similar to those in (i), (ii) above thereafter apply. There 
is also a neat, exact, local solution to the full system (5.2) - (5.4), described in Walton 
(1990), and a possible link with the experimental studies of Holden (1985), where pronounced 
rope-like vortex motions were observed at the edges of hypersonic boundary layers in the 
Mach-number range 11-13. 

The second of the large- interaction structures is rather simpler, its scales being 
implied by the O(M^) thick boundary-layer form. Thus here 

(y,^,T,v,w,p 2 ) ~ AC. (5.5a) 

u = 0(l (5.56) 

leaving the complete nonlinear interaction system (4.2) - (4.4) intact. The only new feature 
here is that the boundary-layer’s normal extent becomes finite, 0 < y < F say, in normalized 
terms. The linearized version then matches up with the linear acoustic- modes analyses of 
Hall and Cowley and Brown and Smith, after which comments as in (i), (ii) again apply. It 
may be significant however that there are infinitely many acoustic modes, and hence possible 
bifurcations, available usually as the flow proceeds downstream, in contrast with the single 
vorticity mode; and the modes present with large wavenumbers, as well as mode-crossing, 
may shed extra light on the nonlinear interaction process. 

(v) Extreme surface conditions, values of the parameters present such as a, C, or pressure 
gradients, could also provide extra insight. 

(vi) Besides the many applications within boundary-layer transition summarized above, 
and the corresponding computational tasks, and nonlinear similarity forms for instance, there 
are also many other flow configurations to which vortex/Rayleigh-wave interaction applies 
in principle (see Section 6 and Appendix B). 

(vii) Finally here, we should mention the possibility that the Smith and Walton (1988, 
1989) breakdown summarized in (3.2) also describes the ultimate behavior of the present 
vortex/ compressible- Rayleigh nonlinear systems, at least for wide vortices. If so, the bound- 
ary layer again separates, effectively, splitting into two increasingly far-apart layers, with 
only slow motion in-between, at a finite distance downstream. That opens up an intriguing 
prospect, namely that a vortex/wave interaction (either as in Section 2 or as in Section 4) 
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can induce separation (as above), which then introduces an extra (inflexional) inviscid mode 
and thence an extra vortex/wave nonlinear interaction, which leads to another separation, 
hence another interaction, and so on. This self-generating process causes the whole flow 
structure locally to cascade into smaller length scales, and thence into substructures. 

6. FURTHER COMMENTS 

The present work extends our previous studies of vortex/wave interaction (Hall and 
Smith (1987-1989)) to the strongly nonlinear regime in which the entire mean-flow profile, 
at any station x, is altered substantially from its undisturbed laminar form. This is for the 
compressible boundary-layer setting in the nonlinear vortex/Rayleigh interaction of Section 
4 (the incompressible version is noted in Appendix B) and for the incompressible case in the 
nonlinear vortex/TS interaction of Section 2, with the corresponding compressible vortex/TS 
interaction following readily from a combination of those two sections (see also Smith and 
Walton (1988, 1989)) and below. The solution properties and suggestions given in Section 3, 
Section 5 hint fairly strongly, we believe, at the potential power of these nonlinear interactions 
in terms of full transition of the flow. In addition, it appears that numerous waves can 
be triggered off (see Section 5), all interacting nonlinearly with the unknown mean flow, 
especially if separation for instance is approached. These and other features of the nonlinear 
interactions found offer exciting prospects for more complete theoretical understanding of 
fully fledged boundary-layer transition. 

Some numerical work for fully nonlinear interaction is described in Section 3, for the 
vortex/TS case, but further concerted efforts on full computations, for both the vortex/TS 
and the vortex/Rayleigh cases, are undoubtedly necessary and these represent a major chal- 
lenge. They should enable quantitative comparisons with experiments and direct numerical 
simulations to be made eventually (an encouraging point being that qualitatively the flow 
structures in Section 2, Section 4 seem to be in line with the numerical-simulation experi- 
ence, e.g., of Kleiser, that many more span wise than streamwise wavenumbers are required to 
accurately describe fully fledged transition, except in its later stages). The pure spatial prob- 
lem, the pure temporal problem, and the combined problem featuring the operator d t + ud x 
are all of interest here since they are likely to provoke different ultimate behaviors and hence 
possibly different views of such phenomena as lambda- vortex formations (see above) and the 
successive collapses in scales referred to in Section 5. The latter in turn may lead on to the 
3D Euler stage locally, cf. Smith (1988) for nonlinear TS transition, allowing comparisons 
with experiments and direct numerical simulations (e.g., Zang, Erlebacher, Hussaini, Kleiser, 
and Biringen) on the later stages of transition, including spikes, intermittency and streak 
production. In all this, however, the computational tasks set by the nonlinear interaction 
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problems posed in Section 2 and in (4.2) - (4.4) for instance seem to present a main hurdle. 

Certain interesting other aspects of the vortex/wave interactions should also be put on 
record here. These are: various limiting or simplified cases, such as those in Section 5; 
similarity solutions (c.f., Section 3, Walton (1990)); the analogy with Gortler- vortex devel- 
opment (based on Hall (1982, 1983) and subsequent works); non-equilibrium critical-layer 
effects (in the context of Section 4); high-frequency properties (for Section 2, Section 4); and 
the implications for by-pass transition where the nonlinear interactions do not start from 
linear small-disturbance growth. 

The nonlinear vortex/wave interactions also have numerous potential applications for 
other types of flows, in all of which the essential ideas of Section 2 and/or Section 4 would 
seem to apply: in wakes, channel flows (Bennett et al. (1988)), pipe flows (Walton, 1990), 
plane Couette flow, water motions, free shear layers, separations, vortex breakdown, where a 
swirl-velocity component is added to the vortex motion (and separation such as in (3.2) could 
correspond to abrupt vortex thickening), and flows over surface roughnesses, for example. See 
also Appendix B. The ideas apply further to length scales other than those taken in Section 
2, Section 4. Again, there may be extensions of interest along the following lines: upper- 
branch flow structures, connected with the high-frequency limit; non-Chapman fluids and 
real-gas effects; surface-cooling effects; external-shock interactions in hypersonic flow, where 
the acoustic modes (Section 5) could play an important part; external pressure-gradient in- 
fluences; and alternate compressible interactions. Nevertheless, the major challenge pressing 
seems to us to be the computational one of accurately solving the vortex/wave nonlinear 
interaction equations set up in Section 2, Section 4, given the encouraging guidelines on the 
strongly nonlinear effects possible. 
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APPENDIX A. THE CRITICAL-LAYER BEHAVIOR 


The critical layer occurs when y = f(x, z)+Re * Y, with Y of 0(1), and the flow properties 
produced are predominantly linear. The solution takes the form 

u = c + Re-*A 1 (a,*)r + . . . + R e~5 A^ (1) + . . . + R e~*(C\U + £ x C/ (2) ) + . . . , (Al) 

« = ... + Re-^A + R e~*(C 2 JjW + C X V {2) ) + . . . + Re"» + . . . , (A2) 

w = ...+ Re^APT* 1 ) + . . . + Re'5(£ 2 F + C X W&) + . . . , (A3) 

p = ... + Re-^jPW + Re"f AP (2) + . . . , (A4) 

where A x is the streamwise vortex shear at the critical layer, the main wave part (superscript 
(1)) depends on the fast scales (X, t) = Re5(x,£) with wavespeed c as in Section 4, and 
the main vortex part (overbarred) is independent of the fast scales. Here we discuss the 
incompressible case with zero pressure gradient, for the sake of clarity; the compressible 
version follows along similar lines. 

The successive balances in mass conservation resulting from (Al) - (A4) give 

= A# (A5) 

uP + vi 2) + wi 1] - fcW™ = 0, (A6) 

-fxX x + V Y + W T = 0, (A7) 

and the x-momentum balances of concern are 

V (1) = (AS) 

A 1 YU i x ) + A x V< a > + W W \ X >Y - W {2) f,\ x = -Pjp + A U$}, (A9) 

< UWUP + + W^ui l) - WWfcUjp > +V\i + WYXu = ... + A V YY . (A10) 

Likewise, the z-momentum balances here become 

XiYW-jp = -Pi 1} + f T P? + A (All) 

< + wWwjp - WWfgWP >=... + AW yy , (A12) 

while the y-momentum equation implies that P^ is independent of Y, and hence equal to 
p^\ and 

P$r 2) = -A x Y’v£ 1) + AVy Y . (A13) 

In (A10), (A12), <> refers to the vortex components only, in the enclosed terms. 
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The wave part in this non-flat critical layer may be analyzed by putting V ^ — fjW^ = 
V a , say, which leaves U^\ V a , satisfying essentially the equations of a linear 3D distur- 
bance in a standard flat critical layer, from (A6), (A9), (All), (A13) and since 
from (A5), (A8). In particular, satisfies a forced Airy equation in the form 

\ 1 YW ( x ) = -Pi 1} A' 1 + A (414) 

which yields the solution 

W W = -(A 2 \ x ia)-^hJ{Y)E + c.c., Y = * Y, (415) 

where J satisfies J" — Y J = 1, J"(±oo) = 0, and can be expressed in terms of the Airy 
function, and E = exp(iaAT — ini). This and the corresponding solutions for V a , U I 1 ' are 
smooth for all Y and satisfy the asymptotic conditions of matching, 

„ Y- 1 {F_ 1 ,G- U H-i\E+c.c.,{V«\wW] -> [Go, H„]E+c.c., as Y - ±oo, 

(-116) 

implied by (4.7), (4.8) (with Appendix B), as required. 

Then the main vortex equation of interest here is that for W, in (A12). Here the behaviour 
in (A16) shows that the nonlinear forcing term on the left-hand side of (A12) decays as Y~ 2 
as Y — + ±oo. So double integration of (A12) with respect to Y produces the response 

W ~ -Btn\Y\ as Y -+ ±oo (A17) 

with B = + (fjHo ~ G Q )H*_\ + (fjHo ~ Go)H-i being found to equal -qA~% as 

defined in (4.4b) or in Appendix B. The logarithmic behavior in (A17) and the associated 
logarithmic terms in U, V at large |K| are responsible for the scale factor C\ in (4.1), as well 
as for the effective slip condition in (4.4b). The condition (4.4b) follows from the match with 
the flow solution in Section 4, since ln\Y\ — > |^nRe + in\y — f\ in effect. 
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APPENDIX B. VORTEX/RAYLEIGH- WAVE INTERACTIONS IN THE IN- 
COMPRESSIBLE CASE 


Many of the extra applications mentioned in the text are concerned more with incom- 
pressible fluids, for which the vortex/wave interaction equations of (4.2) - (4.4) continue to 
hold provided p, /Z are replaced by unity in (4.2a-d), M is replaced by ( u — c) in (4.3a), with 
the Af 2 term omitted, and (4.4b) is replaced by 


<7 = 


A"! 

(uOW 


{}. 


(51) 


where the curly brackets signify the curly-bracketed expression in (4.4b). The governing 
equations in that case are (4.2a-d,h,i), (4.3a-c), (4.4b-d), with the above modifications. In 
addition, the local pressure expansion (4.5a) remains valid provided that in (4.6a-c) Mi,M 2 
are replaced by respectively. Similarly, the local velocity expansions (4.7a-c) remain 

true provided that 


G - 1 = —pi/iaui, H - 1 = ( ftP\ - Poz)/icmi, i = -H^u/ia (52, 3, 4) 


instead of (4.8a-c) in turn; also 


G 0 = -t^G-i/ui - 2p 2 /iau 1 , 


(55) 


Ho = -u 2 H- i/tii + (2fjp2 - Pn)/iaui 


(56) 


in this case. 
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Figure 1 (a) Typical strongly nonlinear vortex/wave interactions in a boundary layer. 

The scales of the slowly varying vortices are shown; the travelling waves present 
are either TS waves (Sections 2, 3) or Rayleigh waves (Sections 4, 5). 
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mini § 




Figure 1 (b) Cross-section 0(Re ^xRe *) in the vortex/ Rayleigh- wave case of Sections 

2 

4, 5, including the critical-layer 0(Re _s ) thick) effect at ^ = /. 


29 












Report Documentation Page 

bc.ce •Vr n - n -sr'3ior 

1. Report No 

NASA CR- 181963 
ICASE Report No. 89-82 

2 . Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

ON STRONGLY NONLINEAR VORTEX/WAVE INTERACTIONS IN 
BOUNDARY-LAYER TRANSITION 

5. Report Date 

December 1989 

6. Performing Organization Code 

7. Authorise 

Philip Hall 
Frank T. Smith 

8. Performing Organization Report No. 

89-82 

10. Work Unit No. 

505-90-21-01 

9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665-5225 

11. Contract or Grant No. 

NAS1-18605 

13. Type of Report and Period Covered 

Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Langley Technical Monitor: 
Richard W. Barnwell 

Final Report 


To be submitted to Journal 
of Fluid Mechanics 


16. Abstract 

The interactions between longitudinal vortices and accompanying waves con- 
sidered hfire are strongly nonlinear, in the sense that the mean-flow profile 
throughout the boundary layer is completely altered from its original undisturbed 
state. Nonlinear interactions between vortex flow and Tollmien-Schlichting 
waves are addressed first, and some analytical and computational properties are 
described. These include the possibility in the spatial-development case of a 
finite-distance break-up, inducing a singularity in the displacement thickness. 
Second, vortex/Rayleigh-wave nonlinear interactions are considered for the com- 
pressible boundary-layer, along with certain special cases of interest and some 
possible solution properties. Both types, vortex/Tollmien-Schlichting and vort- 
ex/Rayleigh, are short-scale/long-scale interactions and they have potential 
applications to many flows at high Reynolds numbers. Their strongly nonlinear 
nature is believed to make them very relevant to fully fledged transition to 
turbulence. 


17. Key Words (Suggested by Author(s)) 

vortex/wave interaction, 
Tollmien-Schlichting 


18. Distribution Statement 

02 - Aerodynamics 
Unclassified - Unlimited 


19. Security Classif. (of this report) 

— * t. ... 

20 Security Classif. (of this page! 

21 No. of pages 

22. Price 

Unclassified 

Unclassified 

36 

A0 3 


NASA FORM 1626 OCT 86 


NAS A -Langley, 1990 




